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I. INTRODUCTION 

The physics of many-particle systems differs fundamentally from the one of a few particles and gives rise 
to new interesting phenomena, such as phase transitions [1, 2] or quantum computing [3-6]. Entanglement 
theory, in particular, appears to have a much more complex and richer structure in the N-partite case than 
it has in the bipartite setting. This is reflected by the fact that multipartite entanglement is a very active 
field of research that has led to important insights into our understanding of many-particle physics [7-14]. 
In view of this, it seems worthy to investigate also the behavior of entanglement measures for large-scale 
systems. Despite the fact that the number of entanglement parameters scales exponentially in the number of 
particles [15], it is sometimes possible to capture the most relevant physical properties by describing these 
systems in terms of very few parameters. 

Recently a duality between highly entangled W states and product states has been established [16]. The 
important class of W states [17] represents a particular interesting set of quantum states associated with 
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high robustness against particle loss and nonlocal properties of genuine entangled multipartite states [18- 
21]. And different experimentally accessible schemes to generate multipartite W states have been proposed 
and put into practice over the years [22-25] 

The duality specifies a single-valued function r of entanglement parameters. We shall refer to r as the 
entanglement diameter, as it will play a crucial role throughout this article. Another reason for the term 
entanglement diameter is that r can be interpreted geometrically as a diameter of a circumscribing sphere. 
The geometrical interpretation and its illustration will be presented in the appendix and now we focus on 
the physical significance of r. 

The entanglement diameter uniquely defines the maximal product overlap and nearest product state [7, 
26-28] of a given highly entangled W state. It has two exceptional points in the parameter space of W 
states. At the second exceptional point the reduced density operator of a some qubit is a constant multiple 
of the unit operator and then the entanglement diameter becomes infinite. The maximal product overlap g of 
these states is a constant regardless how many qubits are involved and what are the values of the remaining 
entanglement parameters. These states are known as shared quantum states and can be used as quantum 
channels for the perfect teleportation and dense coding. Thus the shared quantum states are uniquely defined 
as the states whose entanglement diameter is infinite. 

Furthermore, highly entangled W states have two different entangled regions: the symmetric and asym- 
metric entangled regions. In the computational basis these regions can be defined as follows. If a W state 
is in the symmetric region, then the entanglement diameter is a fully symmetric function on the state pa- 
rameters. Conversely, if a W state is in the asymmetric region, then there is a coefficient c such that the c 
dependence of the entanglement diameter differs dramatically from the dependencies of the remaining co- 
efficients. Hence the point of intersection of the symmetric and asymmetric regions is the first exceptional 
point. It depends on state parameters and its role has not been revealed so far. One thing was clear that the 
first exceptional point does not play an important role for three- and four-qubit W states [29, 30]. 

In this article we show that the first exceptional point is important for large-scale W states. It approaches 
to a fixed point when number of qubits N increases and becomes state-independent(up to 1/N corrections) 
when N S> 1. As a consequence the entanglement diameter, as well as the maximal product overlap, 
becomes state-independent too and therefore many-qubit W states have two state-independent exceptional 
points. The underlying concept is that states whose entanglement parameters differ widely, may neverthe- 
less have the same maximal product overlap and this phenomenon should occur at two fixed points. This 
is an analog of the universality of dynamical systems at critical points. It is an intriguing fact that systems 
with quite different microscopic parameters may behave equivalently at criticality. Fortunately, the renor- 
malization group provides an explanation for the emergence of universality in critical systems [1, 2, 31]. 
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The developed concept distinguishes three classes of W states. The first class consists of highly 
entangled W states which are below both exceptional points and then r varies from r m - m = 1/2 to 
?"o ~ l/\/3 + 0(1/N). We will show that these states are in the symmetric region and their entangle- 
ment diameter is a slowly oscillating function on entanglement parameters. Accordingly, the maximal 
product overlap is an almost everywhere constant close to its greatest lower bound. Similar results have 
been obtained in Ref.[32], where it is shown that almost all multipartite pure states with sufficiently large 
number of parties are nearly maximally entangled with respect to the geometric measure [7] and relative 
entropy of entanglement [8]. We will not analyze rigorously these states since they are too entangled to be 
useful in quantum information theory [33]. 

The second and most interesting class consists of highly entangled W states which are between two 
exceptional points and then r varies from tq to infinity. These states are in the asymmetric region and 
the behavior of the entanglement diameter is curious. We will show that r is a one-variable function in 
this case and depends only on the Bloch vector b of a single qubit. As a consequence g depends only on 
the same Bloch vector too and its behavior is universal. That is, regardless how many many qubits are 
involved and what are the remaining N — 1 entanglement parameters the function g(b) is common. We will 
compute analytically g(b) and thereby find the Groverian and geometric entanglement measures [7, 28] for 
the large-scale W states even if neither the number of particles nor the most of state parameters are known. 

The third class consists of slightly entangled W states which are above both exceptional points. In this 
case the maximal product overlap takes the value of the largest coefficient and these states do not posses an 
entanglement diameter. We will not analyze this trivial case, but will combine the functions g(b) for slightly 
entangled and highly entangled asymmetric W states and obtain an interpolating function g(b) valid for both 
cases. It is in a perfect agreement with numerical solutions and quantifies the many-qubit entanglement in 
high accuracy(Ag/g ~ 1CT 3 at N ~ 10). 

The importance of the interpolating formula in quantum information is threefold. First, it connects 
two quantities, namely the Bloch vector and maximal product overlap, that can be easily estimated in 
experiments [34, 35]. Second, it is an example of how do we compute entanglement of a quantum state 
with many unknowns. Third, if the Bloch vector varies within the allowable domain then maximal product 
overlap ranges from its lower to its upper bounds. Then one can prepare the W state with the given maximal 
product overlap, say g$, bringing into the position the Bloch vector, say g(bo) = go. 

This article is organized as follows. In Sec.II, we review the main results of Ref.[16]. In Sec. Ill, we 
consider two- and three-parameter W states in the symmetric region and show that all of these states are 
almost maximally entangled. In Sec. IV, we consider three- and four-parameter W states in the asymmetric 
region and compute explicitly their maximal product overlap. In Sec.V, we generalize the results of Sec.III 
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and Sec.IV to arbitrary many-qubit W states. In Sec. VI, we discuss our results. In the appendix, we provide 
a geometrical interpretation for the entanglement diameter. 

O. MAXIMAL PRODUCT OVERLAP OF W STATES 

In the computational basis N-qubit W states can be written as 

|W n ) = d|100...0) + O2IOIO...O) + • • • + cjv|00...01), (1) 

where the labels within kets refer to qubits 1,2,...,N in that order. The phases of the coefficients can be 
absorbed in the definitions of the local states = 1,2,..., N) and without loss of generality we consider 
only the case of positive parameters. For the simplicity we assume that is the maximal coefficient, that 
is, cat = max(ci, c 2 , • • • , cat). 

The maximal product overlap g(ip) of a pure state is given by 

g(ip) = max | (i/j\u 1 u 2 ...u n ) |, (2) 

U1,«2,...,«JV 

where the maximization runs over all product states. The larger g is, the less entangled is Hence for a 
quantum multipartite system the geometric entanglement measure E g is defined as 

E g = -2\og g{^). 

The maximal product overlap demarcates three different entangled regions in the parameter space of W 
states: 

1. The symmetric region of highly entangled W states, where g{c\, C2, cat) is a symmetric function 
on all coefficients a. 

2. The asymmetric region of highly entangled W states, where the invariance of g{c\ , C2, cjv) under 
the permutations of coefficients a ceases to be true. 

3. The region of slightly entangled W states, where the inequity g 2 (c\, C2, cat) > 1/2 holds. 

The appearance of the three entangled regions is the consequence of the existence of the two critical 
values for the largest coefficient cat. The first critical value n(ci, c 2 , cat_i) is the solution of 

v / r 2_^ + v / r 2_^ + ... + v / r 2_ t J_ i = (J v-2)r 1> (3) 

which always exists and is unique. Note that the first critical value r\ for the coefficient cat depends on 
the remaining coefficients Cj, i = 1, 2, N — 1 but does not depend on cat. Nonetheless we will use the 
abbreviation ri(cjv) = c 2 , cjv-i) whenever no confusion occurs. 
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The second critical value r$(ci, C2, cjv-i) is given by 

rl = c? + c| + --- + 4_ 1 . (4) 

In what follows we will use the abbreviation ^(cat) = r2(ci, C2, cjv-i) for the simplicity. 

The second critical value is always greater than the first one and thus there are three cases. The first case 
is cat < r\ and the maximal product overlap is expressed via the fully symmetric entanglement diameter 
r(ci, C2, cjv), which is the unique solution of 



- cf + \/r 2 - c 2 + • • • + \/r 2 -c 2 N = (N-2) r. (5) 

Then g is given by 



5 2 = ^v^ ll + A/1-31 [l + X/l-^l-ll + A/l-^l (6) 




and is a bounded function satisfying the inequalities c 2 N < g 2 (c\, C2, cat) < 1/2. 

The second case is r\ < < r^. In this case the entanglement diameter r(c±,C2, cat) is the unique 
solution of 



r2 ~ c i + V r<2 - c 2 + \Jr 2 -c 2 N = (N -2)r (7) 

where only the last radical has the — sign. Then g takes the form 

where again the negative root is taken from the last radical. The expression (8) also has an upper and lower 
bounds and the inequalities c 2 N < g 2 (c\,C2, cat) < 1/2 hold everywhere in the asymmetric region. 
The third case is cat > r2 and g takes the value of the largest coefficient in this case 

g 2 = c 2 N . (9) 

Now g is bounded below and satisfies the inequality g 2 > 1/2. 

Despite the fact that there exist three different expressions for the maximal product overlap it is a con- 
tinuous function on state parameters. Indeed, at cn = r\ both Eqs. (5) and (7) have the same solution 
r = t\ = cat and expressions (6) and (8) for g coincide. At cat — > the solution of (7) goes to infinity, 
r — > oo, and (8) asymptotically comes to (9). At this limit g 2 = c 2 N = r\ = 1/2 and thus the surface 
<? 2 (ci, C2, cat) = 1/2 separates out slightly and highly entangled W states. 
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III. SYMMETRIC ENTANGLEMENT REGION 

In this section we analyze the maximal product overlap of two- and three-parameter W states that belong 
to the symmetric region of entanglement and show that if all coefficients are small, then r is a slowly 
oscillating function close to 1/2. 

A. Two parameter W states 

Equations (5) and (7) are solvable for N = 3 and the answer is [29] 

2R, if c§ < cf + 4 



(10) 



c 3 , if c\>c\ + c\ 



where R is the circumradius of the triangle c\, C2, C3. 

When N > 4 Eqs. (5) and (7) cannot be explicitly solved to give analytic expressions for r in terms of 
the coefficients Ck unless the state posses a symmetry. For example, for N = 4 the equations are solvable if 
any two coefficients coincide and unsolvable if all coefficients are arbitrary [30]. 

However, when N 3> 1 the situation is different. In many cases one can derive approximate solutions 
that quantify the entanglement of W states in high accuracy. We will find such approximate solutions and 
compare them with the exact or numerical solutions. 

Consider first a W states with N = m + k qubits and coefficients 

ci = c 2 = • • • = Cm = a, c m+ i = c m+2 = ■■■ = c m+k = b. (11) 
When m > 1 and n > 1 the state is in the symmetric region and Eq.(5) is reduced to 

mvV 2 - a 2 + kVr 2 - b 2 = (N - 2) r. (12) 

This equation is solvable by radicals. Setting a = cos 9 / \^rn, b = sin e/Vk one obtains 

l)(mcos 2 6> + fcsm 2 6>) + 
16(iV- l)(m- l)(k- 1) 

where 



2 _ 2Nmk - 4(N - l)(mcos 2 6 + fcsin 2 6) + 2mk(N - 2)y/D 



D = 1 - ^— - sin 2 26. (14) 
mk 

At m = 1 or k = 1 the denominator and numerator vanish in Eq.(13), but their ratio gives the correct 
answer. We will not consider this case since it is analyzed in detail in Ref.[30]. 
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If m, k S> 1, then r is almost constant since 



(15) 



The question is when (15) achieves the required accuracy. It can be understood by reference to Fig.l, where 
the 9 dependence of the exact solution (13) is plotted. The graphics show that Ar/r ~ 10 -2 at N ~ 10. 




FIG. 1: (Color online) The plots of the 6 dependence of the exact solution r{ff) for the state (11). The top, middle and 
bottom lines represent the cases (m = 10, k = 10), (m = 12, k — 18) and (m — 30, k = 30), respectively. 
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FIG. 2: (Color online) The maximal product overlap function g 2 (9) at different values of m and k. The axes origin is 
put at the point (0, 1/e) to make it easer the comparison of the exact ant approximate solutions. The top, middle and 
bottom lines correspond to the values (m = 10, k = 10), (m = 12, k = 18) and (m = 30, k = 30), respectively. 




As a consequence of Eq.(15) g 2 is also almost constant and close to its lower bound 1/e [36]. Indeed, 
using approximations 



one obtains 

/ = I + (I)+o(iY (17) 

e V m I \k I 
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The behavior of the maximal product overlap g(9) given by Eqs. (6) and (13) is plotted in Fig.2, which 
shows that Ag/g ~ 10~ 2 at m, k ~ 10. It is difficult if not impossible to observe such small deviations of 
the maximal product overlap in experiments and therefore approximate formulas (15) and (17) have a good 
accuracy when N > 20. 

B. Three parameter W states 

Consider now a three-parameter W state with N = m + k + I qubits and coefficients 

ci = • • ■ = Cm = a, c m+ i = • • • = c m+ k = b, c m+ k+i = • • • = c m+ k+i = c. (18) 

We will analyze the case m,k,l S> 1. Then Eq. (5) can be rewritten as 

m\Jr 2 - a 2 + k\J r 2 — b 2 + l\J r 2 — c 2 = (N — 2) r. (19) 

From the normalization condition ma 2 + kb 2 + lc 2 = 1 it follows that a 2 < l/m <C 1 and similarly 
6 2 , c 2 <C 1. On the other hand (19) shows that r ~ 1, and therefore we can expand the radicals in powers of 
a 2 /r 2 , b 2 /r 2 and c 2 /r 2 . Then 

Again we got the same answer for r, which means that for partitions with large number of qubits r depends 
neither on m, k, I nor on a, b, c. More precisely, r depends only on the expression ma 2 + kb 2 + lc 2 = \ip\ 2 , 
which drops out owing to the normalization condition. 

The equation (19) can be solved explicitly, but the resulting half -page answer is impractical and we will 
compare (20) with the numerical solution instead. For this purpose we use the parametrization 

a = sin 9 cos ip/y/m, b = sinOsmip/^/k, c = cos6. 

The behavior of the numerical solution r(8) of Eq.(19) for various values m, k, I and ip is plotted in Fig. 3. 
The graphics show that the approximate solution is in a perfect agreement with the numerical solution for 

iV»l. 

In summary, in the symmetric region of highly entangled W states the maximal product overlap does not 
depend on state parameters when many qubits are involved. Consider a W state, where n\, n2, product 
vectors in the computational basis have coefficients ci, C2, Cj., respectively. Then g does not depend on 
partition numbers rij or amplitudes q and the approximate solution (15) with the maximal product overlap 
(17) quantifies the entanglement in high accuracy. For example, at N ~ 10 the accuracy is Ag/g ~ 10~ 2 . 
This approximation is true unless the condition rij 3> l(i = 1, 2, k) is violated. What is happening if 
this condition is violated, is analyzed in the next section. 
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FIG. 3: (Color online) The curves show the dependence of the function r(0). The upper, middle and bottom curves 
represent the cases (m = k = I = 10, ip = 7r/4), (m = k = I = 20, ip = 5ir/l2) and (m = 10, k = 20, 1 = 30, ip = 
7r/6), respectively. 



IV. ASYMMETRIC REGION OF ENTANGLEMENT 



In this section we consider three- and four-parameter W states in the asymmetric region and show that if 
one of coefficients exceeds the first critical value r\ , then r is a rapidly increasing function and ranges from 
one-third to infinity when the maximal coefficient ranges from the first critical value to the second critical 
value. 



A. Three-parameter W states 

Consider now the case when I = 1 in (18) 

ci = • • • = c m = a, c m+ i = • • • = c m+ k = b, c m+ k+i = c. (21) 

If c 1, then c/r is small and r is almost constant. This case is analyzed in the previous section and 
now we focus on the case when c/r cannot be neglected. Then either c < r\ or r\ < c < r2- 
When c < n Eq.(5) takes the form 



m 



V^ 2 - a 2 + k\Jr 2 - b 2 + v ' r 2 - c 2 = (N - 2) r. (22) 



The ratios a/r and b/r are small since m, k 3> 1. Hence we expand the radicals in powers of these ratios 
up to quadratic terms and solve the resulting equation. The answer is 



1 !-C 2 / C 2 1-3C 2 . n , 2n o 1 

r= VT3l?' V 1 "^ = T^ 2 "' m ax(a,&)<c<-. (23) 
It is reasonable that r — > 1/2 at c — > 0. 
When c > T\ Eq.(7) takes the form 

m\/ r 2 — a 2 + ksj r 2 — b 2 — V V 2 — c 2 = (N — 2) r. (24) 
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Its approximate solution is 



1 1 



1 



3c 2 



1 



1 



1 2 1 

3 <C <2- 



2 VI -2c 2 ' 
As one would expect, r — > oo at c 2 — > 1/2. 

Surprisingly, both solutions (23) and (25) can be unified to a single solution as follows 

11 C 2 / 2 7 2 \ 2 1 

r = , max a , o < c < -. 

2 VI - 2c 2 V 7 2 



(25) 



(26) 



The question at issue is when (26) gives a required accuracy in the asymmetric region r\ < c < r 2 . 
We compare it with the numerical solutions of (22) and (24) for the values (m = 8,k = 10, a/b = 
0.8, r\ « 0.34) in Fig.4, where the solid line is the plot of (26) and the dashed line is the numerical 
solution. Remarkably, the approximate solution is in a perfect agreement with the numerical one in the 
asymmetric region. 
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FIG. 4: (Color online) Graphic illustrations of the function r(c) for the three- and four-parameter W states. The solid 
curve is the approximate solution (26). The dashed curve is the joined numerical solution of Eqs. (22) and (24). All 
remaining coefficients are well away from the first critical value (ps 0.58) when c varies within the range of definition 
in this case. Accordingly, the state is in the symmetric region when < c < 0.58 and in the asymmetric region 
when 0.58 < c < 0.707. The dotted line is the numerical solution for the state (27). Now another coefficient may 
exceed the first critical value. Therefore there are two first critical values, for the last and the preceding coefficients, 
respectively. The first critical value for the next to last coefficient c is w 0.606 and for the last coefficient d is m 0.59 
which is attained at c = 0.45657. Thus the state is in the symmetric region when 0.45657 < c < 0.606 and in the 
asymmetric region otherwise. Remarkably, the three curves coincide when c > 0.606. 
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B. Four-parameter W states 

However, there are W state that are outside the realm of the model sketched in the previous subsection. 
These are states with few (at most three) coefficients close to the first critical value n ~ 1/ v3- In this case 
these coefficients are not small and the resulting r should has a different behavior. 

Notice, two coefficients cannot exceed the first critical value simultaneously. But we can construct W 
states whose coefficients depend on a free parameter in such a way that at one value of the free parameter 
the last coefficient exceeds the first critical value and at another value of the free parameter the preceding 
coefficient exceeds the first critical value. Below we construct an illustrative example of a such state and 
analyze its entanglement diameter. 

An example is the 19-qubit four-parameter W state with coefficients 

Cl = . . . = Cy = a , c 8 = • ■ ■ = ci 7 = b, c 18 = c, ci 9 = d. (27) 
For the normalized states we can use free parameters ip, k and c as follows 

2 COS 2 Lp 2 2 sm2 9?/-, 2\ r2 k ~ 1 ^ 2^ 

= 7k c )• b = -ior (1 - c >' d = — {1 - c ) - 

Now we analyze the function r(c) at k = 1.8, ip = 7r/4. 

1. The next to last coefficient c coincides with its first critical value r\{c) at c pa 0.606 , that is, the 
solution of the system 



7y r\ — a 2 + 10y r\ — b 2 + J r\ — d 2 = \7r\ and r\ = c 

is ri = c pa 0.606. Then r(c) should range from ri(c) to infinity when c ranges from r\{c) to 1/2 
and should has a vertical asymptote at c 2 — > 1/2. 

2. The last coefficient d coincides with its first critical value r\(d) at d pa 0.593, that is, the solution of 
the system 



7y r\ — a 2 + 10y r\ — b 2 + W r\ — c 2 = \7r\ and r± = d 

is n = d pa 0.593. Note that at this point c pa 0.45657. Then r should increase when d ranges from 
r± (d) to (i max . But the maximum value of d is less than the second critical value since <i 2 nax = d?{c = 
0) = (k — l)/k = 4/9 < 1/2. Therefore r should be bounded above in the interval [r\(d), d max ] 
and attain a maximum at d max . As d is a decreasing function on c, r should attain a maximum at 
c = and then decrease when c ranges from to 0.45657. 
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3. The state is in the symmetric region when d < r±(d) and c < ri(c). Hence r(c) should be minimal 
and nearly constant when 0.45657 < c < 0.606. 

The dotted line in Fig.4 represents the c dependence of the function r(c). It agrees completely with the 
above analyze. 

The main point is that all the three curves coincide when c > ri(c). In the next section we will show 
that this is not accidental and the curves must coincide. In this context the equation (26) is a surprising 
result. The quantity r, as well as the maximal product overlap g, depends from c only. The rest of the state 
parameters appear in (26) in the combination \ip\ 2 — c 2 and drop out by the normalization condition! 




FIG. 5: (Color online) The plots of the function g(c). The solid line is the approximate solution (28), the dashed 
and dotted lines are the numerical solutions for the states (21) and (27), respectively. The curves may have different 
behaviors when cjv < T\, but coincide when cjv > r%. 



Furthermore, we can derive an analytic expression for the maximal overlap. Using approximations (16) 
one obtains 

5 2 (c) = (1 _ c 2 )e -(l-2^)/(l-^)_ (2g) 

The behavior of the function g(c) is shown in Fig.5. The solid line is the curve (28), the dashed curve is 
the numerical solution for the state (21) and the dotted line is the numerical computation for the state (27). 
They all coincide when c > r\(c). 

For highly entangled states the maximal product overlap ranges from its lower to the upper bound when 
c ranges from r\ to r2- On the other hand the Bloch vector b of Nth qubit is collinear with axis z and 
b z = 1 - 2c 2 . Thus g is a one-variable function on b z and one can vary the entanglement of the multiqubit 
W state by altering the Bloch vector of a single qubit. The remaining qubits should be present in order to 
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create an entanglement, but their individual characteristics do not play any role within the domain — 1 < 
b z < 1 — 2r 2 , N S> 1. These qubits are just spectators, they should appear in the W state, but have no 
influence on the entanglement of the state. 

V. GENERAL CASE 

The results of the previous sections are based on the fact that the entanglement diameter r is bounded 
below. In the symmetric region it is rigidly bound by the following theorem. 
Theorem 1. If r is a solution of Eq.(5), then 

\ < r 2 < \ (29) 

Proof. Note that 




By summing over i the above inequality and using (5) and the normalization condition one obtains 

< 2(n - n + 2) = 4. 

Hence r 2 > 1/4. Next, from x < y/x for < x < 1 it follows that 

n ( 2 \ n I ?~ 1 
1=1 v 7 1=1 

that is, r 2 < 1 /2. 

The inequalities (29) allow us to understand the behavior of g of arbitrary N-qubit W states in the 
symmetric region. Indeed, in this region cf ~ 1/N and therefore cf/r 2 <C 1. Then one can expand the 
radicals in (5) and obtain 

N-±~N-2, 

which generalizes (5) and (6) to arbitrary W states with cat <C 1. 

In Eq.(18) we have chosen equal coefficients in order to reduce the number of independent parameters 
and make it easier the analyze. Now Theorem 1 states that it is irrelevant whether some coefficients coin- 
cide. Decisive factor is that the coefficients Cj are small(~ 1/y/N). Then the ratios a/r are small since r 
is bounded below (~ 1/2) and we can keep first nonvanishing orders of these ratios. Surprisingly, all these 
ratios are combined in such a way that they yield the Euclidean norm of the state function and the final 
answer becomes independent on the state parameters as well as the number of particles involved. 
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In the asymmetric region the entanglement diameter r should has a lower bound but has not an upper 
bound since r — > oo at C2 — > T2- One may expect that the lower bound of r in the asymmetric region 
coincides with the upper bound of r in the symmetric region. But the following theorem shows that this is 
not the case. 

Theorem 2. If r is a solution of Eq.(7), then 

r 2 > i (30) 

Proof. We use the same technique, namely 

N-l 2 2 N ~ 1 



— = y^ + ^<V2(i 

>» 2 ^^^^^ [y* 2 lyt 2 ^^^^^ 1 



or 




1 / c 2 c 2 

- T <2-2;/l--^ + -^<3 since cyv < r. 

This bound, as well as bounds (29), is tight, for example, r 2 — » 1/3 at c 2 — )• 1/3 in (26). 

Theorem 2 explains why the asymmetric approximate solution (26) fits the numerical date more quickly 
(N ~ 10) than the symmetric one (15)(iV ~ 20). First, the lower bound of r is greater in this case. 
Second, since cat is greater(cAr > r{) the remaining coefficients should be smaller due to the normalization 
condition. These two factors together make the ratio Cj/r smaller. Hence the approximate solution should 
has a better agreement with the exact one. Aside from that, r is a fast increasing function and goes to the 
infinity unlike to the symmetric case. Hence the values of the coefficients q become irrelevant when r>l. 

In fact there is no W state in the asymmetric region that differs markedly from the above model when 
many qubits are involved. The following theorem completes the proof that in the asymmetric region the 
maximal product overlap is a one-variable function. 

Theorem 3. If cjv = r\, then 

^4 + 0(1) (3D 

Proof. Note that on the boundary of the symmetric and asymmetric regions r = r\ = cn and therefore 
ff > 1/3. Expanding the radicals in (3) in powers of cf/rf one obtains 

which gives (31). 

Now we are ready to explain what is happening in the asymmetric region. 
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1. When many qubits (TV S> 1) are involved the first critical value depends neither the number of qubits 
nor the state parameters and is a constant, r\ » 1/ \/3- 

2. Regardless what is happening in the interval < c^r < T\ all functions r(c) must converge to the 
point r (1/ y/3) w 1/ \/3- This is the effect of the first critical value. 

3. All functions r(c) have the the same vertical asymptote, namely, r(c) — > oo at c — > l/v2. This is 
the effect of the second critical value. 

These statements together give no chance to differ markedly exact and approximate solutions in the 
asymmetric region. In conclusion, when N S> 1, everywhere the maximal product overlap of W states is 
governed by the smallest b z among the z components of the Bloch vectors. Using approximations 

1 



2 V 1 + V 1 "^ 
and equations (9) and (26) one obtains 

g\N » 1) = 



e _ c 2 /4r 2 



1,2,- •• ,N-\ 



1+6* r 
2 e 

l-b z 
2 ' 



2.b z . 

if 0<b z <\ 



if 6 7 < 



(32) 



Graphic comparison of the interpolating formula and numerical computation of g is shown in Fig. 6, where 
the b z dependence of g is plotted for N = 10. The solid and dashed lines represent the interpolating function 
(32) and numerical computation, respectively. 

g 1 



v 1.0 




\ 0.9 




X. 0.8 




>v 0.7 




>v 0.6 




0> 




0.4 





FIG. 6: (Color online) The maximal product overlap g as a function of z component of the Bloch vector b z . The solid 
line is the interpolating formula (32). The dashed line is the numerical computation for a 10-qubit W state. 



We did not plotted numerical results for different states because different curves overlap and become 
indistinguishable. We failed to find the states for which the numerical results markedly differ from the 
plotted one provided N 3> 1 holds. 
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VI. DISCUSSION 

The main result of this work is the formula (32). First, it shows that sometimes the characterization 
and manipulation of the entanglement of many qubit states is a simple task, while the case of few or several 
qubits is a complicated problem. Second, it states that when N 3> 1 the maximal product overlap of W states 
is universal in the asymmetric and slightly entangled regions and the only exceptions are W states in the 
symmetric region that are almost maximally entangled states. Then a question arises: Why do the maximal 
product overlaps of the different W states far apart from the exceptional points have the same behavior? 
Perhaps the reason is that these states are all W-class states. Classification of entangled states explains 
that pure states can be probabilistically converted to one another within the same class by stochastic local 
operations and classical communication [17, 37, 38]. And one can assume that large-scale systems within 
the same class have the feature, aside from the interconvertibility, that their entanglement is universal. An 
argument in favor of this assumption is that the geometric measure of entanglement [7], the relative entropy 
of entanglement [8] and the logarithmic global robustness [9] are related by bounding inequalities and, 
moreover, the relative entropy of entanglement is an upper bound to entanglement of distillation. Hence 
it is unlikely that these measures may exhibit contradicting results and each of them predicts its own and 
very different entanglement behavior of large-scale W-states. If this argument is true, then entanglement of 
large-scale states within the same class is universal. However, states from the different classes may exhibit 
different behaviors. By no means it is obvious, and probably not true, that the maximal product overlap of 
GHZ-class states should have a behavior similar to that of W states. 

Another possible explanation is that the universality of the maximal overlap of large scale W states is 
the inherent feature of the geometric entanglement measure rather than the inherent feature of quantum 
states. If it is indeed the case, then a reasonable question is the following: do the exceptional points really 
exist or they are just the fabrication of the geometric entanglement measure? In this context the second 
exceptional point is a fundamental quantity. Indeed, there are states applicable for the perfect teleportation 
and dense coding and these states all should possess the same amount of entanglement. Hence there is 
an specific entanglement point(infinite entanglement diameter in the case of the geometric measure) that 
can be associated with the exceptional point. And one can assume that the second exceptional point is a 
property of quantum states rather than a property of the maximal product overlap. And how about the first 
exceptional point? Unfortunately, we do not know any strong arguments in favor of it. In order to clarify the 
existence or nonexistence of the first exceptional point, as well as the second exceptional point, one has to 
analyze another reliable entanglement measure, say relative entropy of entanglement [39], and see whether 
it possesses exceptional points. 
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Appendix A: Geometrical interpretation of the duality 




FIG. 7: (Color online) The geometrical interpretation of the duality for three-qubit W states. Mutually perpendicular 
bold lines OX, OY and OZ are coordinate axes and OO' is an arbitrary direction. OCx , OCy and OCz are mirror 
images of the line OO' in respect to the three axes. The points Cx, Cy and Cz are intersections of these lines with 
the sphere uniquely defined by the two conditions: its center lies on the line OO' and its diameter OD = r is the sum 
of the lateral sides of the upperpyramid (with the apex O and base CxCyCz)- Now the direction cosines (and sines) 
of the vector are coefficients of the local states \m) in a computational basis. And the lateral sides of the lower 
pyramid (with the apex D and base CxCyCz) are the coefficients of a 3-qubit W-state in the same basis. Thus each 
direction singles out a product state and a W state and thereby establishes a correspondence among them. 



The nearest product state \ui) ® \w) 8> • • • ® \un) of the W state (1) can be parameterized as follows 



7T 



where 



\u k ) = sin6> fc |0) +cos# fc |l), 0< k < -, k = 1, 2, N, 



cos 2 9\ + cos 2 02 + • • • + cos 2 0N = 1- 



(Al) 



(A2) 



Thus the angles cos 9 k define a unit N-dimensional vector in Euclidean space. They satisfy the equalities 

1 _ sin 20\ sin 2^2 sin2#jv 
r Cl C2 cn 



(A3) 
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These equalities can be interpreted as trigonometric relations for the right triangles with hypotenuses r, 
angles 26 k, opposite legs c& and adjacent legs Jr 2 — c 2 ,. If 20 k > ir/2, then one takes the angle it — 29 ^ 
instead. All of these triangles has the same hypotenuse r and therefore can be circumscribed by a single 
sphere with the diameter r. The final picture represents two inscribed N-dimensional pyramids with a 
common base and lateral sides c\, C2, cjv and V 2 — c\ , \Jr 2 — c 2 , \Jr 2 — c 2 N , respectively. The case 
N = 3 is illustrated in Fig.7. 
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